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Abstract. IceTop, the surface component of the IceCube detector, has been used to measure the energy 
spectrum of cosmic ray primaries in the range between 1.58 PeV and 1.26 EeV. It can also be used to study 
the low energy muons in air showers by looking at large distances (> 300nr) from the shower axis. We will 
show the muon lateral distribution function at large lateral distances as measured with IceTop and discuss the 
implications of this measurement. We will also discuss the prospects for low energy muon studies with IceTop. 


It is well known that the muon content of an air shower, 
together with a measure of its electromagnetic component, 
can be used to estimate the energy and mass of its pri¬ 
mary m. The main issue with the use of the muon content 
as an estimate of primary mass is the possible systematic 
differences between simulated and real air showers, aris¬ 
ing from the lack of knowledge of high energy hadronic 
interactions. An excess number of muons has been re¬ 
ported by the HiRes-MIA and the Pierre Auger collabo¬ 
rations El a. Understanding the systematic differences 
in air shower muon content between simulations and data 
is one of the pressing issues in the physics of very high 
energy cosmic rays. 

The IceTop detector is sensitive to the low-energy (E 
> 200 MeV) muon component of air showers. Generally 
speaking, at low zenith angle and close enough to the air 
shower axis, the signal from muons is overwhelmed by the 
signal from the electromagnetic (EM) component of the 
air shower (electrons, positrons and gamma-rays). This 
holds true in the zenith angle and lateral distance ranges 
that have been used in the cosmic ray spectrum determi¬ 
nation with IceTop ( 6 < 40°, r < 220 ■ ( E/PeV )^ 4 m). 
We are interested in determining the average muon Lateral 
Distribution Function (LDF) at lateral distances where the 
signal from muons becomes significant. The lateral dis¬ 
tance of any point is defined as the closest distance from 
the point to the shower axis. 

The estimation of the muon lateral distribution hinges 
on being able to distinguish the signal produced by single 
muons from the signal produced by EM particles. IceTop 
detector stations accomplish this by virtue of their size. 
A typical electron leaves a signal proportional to its track 
length through the detector, which is on the order of 10 
cm. A typical muon crosses the detector, producing a track 
typically longer than 90 cm. After describing the general 
features of IceTop signals in Section [T| we will describe 
a method for measuring the number of muons in Section 
[2] This method does not rely on counting muons, but it is 


a fit to the signal distributions within small bins of lateral 
distance, primary energy and zenith angle. 

1 General Features of IceTop 

The IceCube detector consists of two major components. 
It can measure air showers on the surface with IceTop, 
high energy muon bundles with the in-ice detector, and 
both components in coincidence provided that the air 
shower axis goes through the in-ice detector. In what fol¬ 
lows, we will consider the specific characteristics of Ice¬ 
Top that are relevant for measuring the low-energy muon 
component of air showers. A more detailed description 
of IceCube and IceTop has already been presented else¬ 
where 0. 

IceTop is an air shower array consisting of 81 stations 
forming a triangular grid with a separation of 125 m in 
its completed configuration. The results presented here 
were obtained with data collected between June 1st 2010 
and May 13th 2011, when IceTop consisted of 73 sta¬ 
tions. It is located above the deep IceCube detector at 
the geographical South Pole, covering an area of roughly 
one square kilometer. Each station consists of two ice 
Cherenkov tanks separated by ten meters. Each tank con¬ 
tains two Digital Optical Modules (DOMs) with a 10 inch 
photomultiplier tube (PMT) and electronics for signal pro¬ 
cessing and readout. A discriminator trigger occurs when 
the voltage in one of the DOMs in a tank has passed 
the discriminator threshold. The total charge collected at 
the PMT’s anode, after digitization and baseline subtrac¬ 
tion, constitutes the tank’s signal. The tanks register sig¬ 
nals ranging from 0.2 to 1000 Vertical Equivalent Muons 
(VEM). A Hard Local Coincidence (HLC) occurs when 
there are discriminator triggers in two neighboring tanks 
within a time window of 1 |is. If there is a discriminator 
trigger but not an HLC, the result is a Soft Local Coinci¬ 
dence (SLC). 

All analyses of IceTop data up to now have only con¬ 
sidered signals where both tanks in a station pass the 
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threshold, or HLC signals. In this contribution we also 
consider SLC signals, where the partner tank within the 
station did not have a discriminator trigger. SLC signals 
occur at large lateral distances, where the triggering prob¬ 
ability is smallest. An example of the lateral distribution 
of SLC and HLC signals from experimental data is shown 
in Figure [T] 

The properties of the primary cosmic ray are recon¬ 
structed by fitting the measured signals with a Lateral Dis¬ 
tribution Function (LDF) which includes an attenuation 
factor due to the snow cover on top of each tank. For a 
given primary energy and arrival direction, the observed 
lateral distribution of signals used in this reconstruction is 
shown in Figure[2] The signal times are fitted with a func¬ 
tion describing the shape of the shower front. The primary 
energy is given by the shower size S 125 , defined as the sig¬ 
nal at a lateral distance of 125 m. The resulting cosmic ray 
energy spectrum measured with IceCube/IceTop has been 
presented previously |5}|2|]. In the shower reconstruction 
just outlined, only HLC signals were considered. 

The distinction between SLC and HLC signals pro¬ 
vides a natural way to identify tanks where we expect 
to see a significant muon contribution. Generally speak¬ 
ing, we expect that signals at large lateral distances will 
be mostly due to muons, whereas the signals at short 
lateral distances will be mostly due to electrons and y- 
rays. Instead of relying on simulations for selecting sig¬ 
nals where the contribution from muons dominate, we will 
select tanks where SLC signals are most likely. Specifi¬ 
cally, we will restrict ourselves to tanks beyond a lateral 
distance at which SLC signals amount to 50% of the total 
number of signals. The exact lateral distance at which this 
happens depends on the energy of the primary. The depen¬ 
dence on zenith angle is very weak. Therefore, we chose 
to use the value corresponding to 30°. 

The statistical distribution of signals from EM parti¬ 
cles will roughly mimic their energy distribution, with a 
mean signal that corresponds to a few tens of centime¬ 
ters of track length inside the tank. On the other hand, 
the distribution of signals from muons is mainly deter¬ 
mined by the geometry of the tank. The signal distribu¬ 
tions produced by single muons are obtained using the 
Geant4 toolkit ( 8 ). Example distributions at various in¬ 
cident angles are displayed in Figure [3] The distribution 
is clearly not symmetric. The peak of the distribution cor¬ 
responds to muons that enter through the top of the tank 
and exit through the bottom. By definition, the peak po¬ 
sition for vertically through-going muons is one Vertical 
Equivalent Muon (VEM). For muons arriving at a zenith 
angle 6, the peak is at 1 / cosW). The flat part at low signals 
corresponds to muons with a short track through the tank, 
what we call corner clipping muons. At large angles, few 
muons go through top and bottom. For an integer number 
of muons, the signal distribution is just the multiple auto¬ 
convolution of the single-particle distribution. An exam¬ 
ple of this is displayed in Figure [4] The total signal dis¬ 
tribution for an expected average number of muons (N^) 
is given by a linear combination of the signal distributions 
for integer numbers of muons: 



Figure 1. Lateral distribution of tanks in showers arriving with 
zenith angles less than 6° and with energies between lOPeV and 
12.6 PeV. The tanks are classified depending on whether they 
registered SLC or HLC signals (described in the text). The upper 
curve includes all tanks, including tanks that registered no signal 
during the event. In determining the muon content, we will only 
consider tanks at lateral distances larger than where the SLC and 
HLC lines cross. 
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Figure 2. HLC signal distribution as a function of lateral dis¬ 
tance for air showers with energies between 4 and 5 PeV, and 
zenith angle between 28° and 32°. This distribution is a 2d his¬ 
togram that includes all HLC signals from all events with the 
given energy and arrival direction. 
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2 Determining the Muon Lateral 
Distribution 

An example of the observed distribution of IceTop signals, 
corresponding to air showers with fixed energy and zenith 
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Figure 3. Signal distributions resulting from simulating the de¬ 
tector response for single muons arriving at fixed zenith angles 
from 0° to 57°. 


Figure 5. Signal distribution as a function of lateral distance for 
air showers with energies between 4 and 5 PeV, and zenith angle 
between 28° and 32°. This distribution includes all signals (SLC 
and HLC). The vertical lines mark the distances at which the 1-d 
histograms in Figure[6]were made. 



Figure 4. Signal distributions resulting from simulating the de¬ 
tector response for an integer number of muons arriving arriving 


at 10°. 


angle, can be seen in Figure [5] where we now included 
the SLC signals. At large distances, there are two dis¬ 
tinct populations. One population is the continuation of 
the main distribution at smaller distances, which roughly 
follows a power law, where the electromagnetic compo¬ 
nent of the shower dominates. The other population, with 
signals around 1 VEM, is made up mostly of tanks hit by 
one or more muons. These two populations are clearly 
seen in Figure [6] where we show the histograms of signals 
registered at selected fixed lateral distances. 

The first population corresponds to tanks hit by no 
muons. The entire signal in these tanks is produced by 
electrons or y-rays. We approximate this distribution by 
a power-law multiplied by a function that describes the 
trigger probability. The trigger probability can be de¬ 
scribed by a sigmoidal function of the logarithm of the 
charge, centered at 0.25 VEM and with a width of 0.14. 


r = 297 m 



r = 363 m 



Figure 6. Vertical slices of the 2-d histogram in Figure [5] Each 
histogram corresponds to a vertical line in Figure[5] The fit lines 
are not shown at 297 m because at this point the assumption made 
of the shape of the EM signal, that the peak EM signal is below 
threshold, is not valid. 


Clearly, this approximation only works well at large lat¬ 
eral distances, where the mean expected signal is well be¬ 
low threshold and we thus look into the tail of the signal 
distribution. The second population, with a peak around 
1 VEM/cos 6, can be described by the contributions of 
tanks hit by an integer number of muons. The charge 
distribution obtained from simulations was described in 
Section |T] and is given by Equation [T] This distribution 
is smeared and shifted to account for a very small contri¬ 
bution from electrons, positrons, and y-rays. 

We determine the muon lateral distribution function 
by independently fitting the charge distributions at fixed 
energy, zenith, and lateral distance, like the distributions 
shown in Figure [6] The result of the fit is the number of 
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Figure 7. Reconstructed average lateral distribution of muons for 
air showers arriving at zenith angles of less than 6 ° and selected 
S {25 values. 



Figure 8. Reconstructed average lateral distribution of muons 
for air showers arriving at zenith angles between 28° and 31°, 
and selected S 125 values. 


tanks hit by at least one muon. This, together with the total 
number of tanks located at that distance (the upper curve 
in Figure[I]i, provides an estimate of the probability that a 
tank is hit by one or more muons, which leads to the mean 
number of muons (A^): 

Ntanks with muons , -/is/ \ 

Pa hit = — - - = 1 - e ' ( 2 ) 

™ all tanks 

The mean number of muons is divided by the cross- 
sectional area of the tanks to yield the muon density at 
that location. In doing this, we assume that the direction 
of motion of all the muons coincide with the reconstructed 
air shower direction. 

The resulting muon density for each lateral distance 
and some of the energies considered, at two different 
zenith angles (0° and 31° from the vertical), are shown in 
Figures [7] and [ 8 ] At this point we should remember that all 


the points in these figures are chosen so the number of SLC 
signals at any of the considered lateral distances are at least 
50% of the number of signals at that lateral distance. The 
smaller and fainter markers correspond to tanks located at 
lateral distances in which SLC signals amount to less than 
80% of all signals at that distance. Clearly, increasing this 
lateral distance cut enhances the muon’s contribution to 
the signal, but at the same time limits the accuracy of the 
lateral distribution fit. The optimization of this cut remains 
under study. 

Each LDF can be described by the following func¬ 
tion 0 : 


Pn(r) = Pit( r 0 ) 



320 m + r \~ r 
320 m + ro ) 


(3) 


which displays the same functional form as Greisen’s 
function, with the first exponent of r fixed to -3/4. The pa¬ 
rameters that are fitted are y and /j /( (/- 0 ), which represents 
the density of muons at a given lateral distance ro. We 
have chosen 600 meters as the value for ro. The choice of 
ro is arbitrary and is motivated by the fact that the signals 
shown in Figures [7] and [ 8 ] are at lateral distances between 
300 and 1000 meters. It is also convenient because this is 
the lateral distance at which previous experiments, notably 
Akeno and HiRes-MIA j2][T0!!]Ti, reported their results. 
It must be noted that the optimum value for ro depends on 
energy and on the lateral distance cut. The two parame¬ 
ters, p^ro ) and y, potentially depend on energy and zenith 
angle. In the following section we will show how p M (r {) ) 
depends on energy in the case of vertical air showers. 


3 Results and Discussion 

The resulting muon density at 600 m from the shower axis 
is displayed in Figure [9] for the case of vertical events. 
Here one can see that the muon density is in the same or¬ 
der of magnitude as the Hires-MIA experiment J2]. One 
remarkable result is that the dependence on energy agrees 
very well with the Akeno air shower array result mm, 
which reported a power law dependence with an index of 
0.83. 

The interpretation of any apparent differences in the 
absolute scale in Figure[9]requires some care. We must re¬ 
member that the MIA array was not located at the same 
atmospheric depth as IceTop (the depths are 860 and 
680g/cm 2 respectively), which means that the densities 
need to be corrected to account for attenuation in the atmo¬ 
sphere. Additionally, the detection thresholds are slightly 
different. However, it is very encouraging that the apparent 
offset in Figure [9] is on the order of the expected attenua¬ 
tion. 

At this moment we do not believe that any deviation 
from a power law behavior in Figure [9] should be given 
much importance. The cause for such deviations can still 
be due to limitations of our method or the choice of lat¬ 
eral distance cut. These and other systematic effects are 
currently under study. 

A more detailed comparison with the expectation from 
simulations will remain for future contributions, after we 
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Figure 9. The energy dependence of the muon density p^(600) 
for vertical showers. Also shown is the result from the HiRes- 
MIA collaboration m. 


have taken into consideration all relevant systematic ef¬ 
fects. In this contribution we note that the statistical er¬ 
rors in the present measurement are significantly smaller 
than any other previous measurement. We expect to make 
a statement on the absolute scale of the muon density at 
large lateral distances in the near future. 

We are considering improvements to the current anal¬ 
ysis. The main limitation at energies larger than ~ 30PeV 
arises from our selection of signals beyond a given lat¬ 
eral distance cut, together with the finite size of the array. 
These two constraints limit the lever arm in the lateral fit. 
One of the possibilities we are considering is the exten¬ 
sion of this analysis to include air showers that go through 
the detector deep in the ice but not through IceTop, which 
would increase our lateral distance range. We could also 
derive improvements from changes in our lateral distance 
cut. 

An obvious application of this measurement will be 
the improvement of our event-by-event reconstruction to 


include a muon LDF together with an EM LDF, which 
we expect will increase our sensitivity to the mass of the 
primary. This method will provide a measure of primary 
composition that is independent of the current approach of 
measuring the high-energy muon bundles deep in the ice 
(Efj <; 300 GeV). We expect that a comparison of these 
two independent measures can help in identifying system¬ 
atic effects arising from different hadronic models used in 
Monte Carlo simulations of air showers. 
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